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Abstract. The main scientific goal of the ICANOE detector(l) is the one of elucidating in a comprehensive way the 
pattern of neutrino masses and mixings, following the SuperKamiokande results and the observed solar neutrinos deficit. 
To achieve these goals, the experimental method is based upon the complementary and simultaneous detection ofCERN 
beam (CNGS) and cosmic ray (CR) events. For the currently allowed values of the SuperKamiokande results, both 
CNGS and cosmic ray data will give independent measurements and provide a precise determination of the oscillation 
parameters. Since one will observe and unambiguously identify v e , and v x components, the full (3 x 3) mixing matrix 
will be explored. 



INTRODUCTION 

The reference mass for underground detectors is now 
set by the operating SuperKamiokande(2) detector, which 
is of the order of 30 ktons. However the rather coarse na- 
ture of the Cherenkov ring detection is capable to recon- 
struct only part of the features of the events. 

New generation underground experiments are now 
facing new challenges, for which novel and more power- 
ful technologies are required, with respect to the existing 
detectors: 

1. The long baseline accelerator neutrino oscilla- 
tion experiments require, with respect to exist- 
ing short baseline detectors (like NOMAD(3) and 
CHORUS(4)), a large increase of the detector fidu- 
cial mass (about 3 ton in NOMAD), in order to cope 
with the flux attenuation due to the distance. In or- 
der to perform a comprehensive program on neutrino 
oscillations the fiducial mass of the detector must be 
increased to several ktons. Moreover, the detector 
must be able to tag efficiently the interaction of v e 's 
and v x 's out of the bulk of v^, events. This requires 
a detailed event reconstruction that can be achieved 
only by means of a high granularity detector. 

2. Likewise, the comprehensive investigation of atmo- 
spheric neutrinos events, in order to reach the level 
of at least one thousand events/year, also requires 
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a fiducial mass of several ktons. The capability to 
observe all separate processes, electron, muon and 
tau neutrino charged currents (CC) and all neutral 
currents (NC) without detector biases and down to 
kinematical threshold is highly desirable. 

3. Nucleon decay: because of the already very high 
limit on the nucleon lifetime (> 10 32 years in most 
of the decay channels), a modern proton decay 
detector should have an adequately large sensitive 
mass. 

Event imaging should be provided by a modern bubble 
chamber-like technology since (1) it has to be able to pro- 
vide high resolution, unbiased, three dimensional images 
of ionising events; (2) it has to provide an accurate mea- 
surement of the basic kinematical properties of the parti- 
cles of the event, including particle identification. (3) it 
has to accomplish simultaneously the two basic functions 
of target and detector. 

In the ICANOE design a fully sensitive, bubble 
chamber-like detector will permit discovery limits at the 
few events level and a much more powerful background 
rejection. A detector of this kind, already at the level of a 
few ktons of mass will be fully competitive with the po- 
tentialities of SuperKamiokande and in several domains 
will permit to extend much further the investigations. 

The ICANOE detector fruitfully merges the superior 
imaging quality of the ICARUS technology(5) with the 
high resolution full calorimetric containment of NOE(6), 
suitably upgraded to provide also magnetic analysis of 
muons. It has a modular structure of independent super- 
modules and is expandable by the addition of such super- 
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FIGURE 1. Perspective view of the baseline detector with 4 
supermodules. 



modules, each consisting of a low density 1 .9 kton liquid 
target and of a high density 0.8 kton active solid target. 

The superior quality of the event vertex inspection 
and reconstruction of the liquid Argon is ideally com- 
plemented by the addition of the external module capa- 
ble of magnetic analysis of the muons escaping the LAr 
chamber. Bubble chambers have in fact often been very 
similarly complemented in the past by external identi- 
fiers. An iron muon tracking spectrometer would fulfill 
this job, but it would also introduce in between adjacent 
liquid argon volumes a blind region incapable of giving 
information on the energy and on the nature of the es- 
caping particles. A sensitive magnetized calorimeter ap- 
pears therefore as an ideal containment module to be in- 
terleaved between adjacent liquid Argon volumes. 



OUTLINE OF THE ICANOE DETECTOR 

The ICANOE layout (Figure 1) is similar to that of a 
"classical" neutrino detector, segmented into almost inde- 
pendent supermodules. The layout of the apparatus can 
be summarized as follows: 

• the liquid target, with extremely high resolution, 
dedicated to tracking, dE/dx measurements, full 
e.m. calorimetry and hadronic calorimetry, where 
electrons and photons are identified and measured 
with extremely good precision and n/p,K and p sep- 
aration is possible at low momenta; 

• the solid target, with good e.m. and hadronic res- 
olution, dedicated to calorimetry of the jet and a 
magnetic field for measurement of the muon features 
(sign and momentum); 



• The supermodule, obtained joining a liquid and 
solid module, which constitutes the basic module of 
an expandable apparatus. A supermodule behaves 
as a complete building block, capable of identify- 
ing and measuring electrons, photons, muons and 
hadrons produced in the events. The solid, high 
density sector reduces the transverse and longitudi- 
nal size of hadron shower, confining the event (apart 
from the muon) within the supermodule. 

The ICANOE SuperModule is, according to the 
present design, composed by a liquid Argon module, with 
18. Ox (11.3) 2 m 3 of external dimensions and 1.4 kton 
(1.9 kton) of active (total) mass, and a magnetized 
calorimeter module, with 2.6 x 9 2 m 3 of external dimen- 
sions and 0.8 kton of mass. The depth consists of 2 m for 
the calorimetric unit, corresponding to 7.4 A,,-,,, and 59 Xq 
and 0.6 m for the tracking units (the interleaved planes of 
tracking chambers). 

At this stage, four SuperModules with a total length 
of the experiment of 82.5 m and a total active mass of 
9.3 kton fully instrumented are being considered for the 
baseline option. 



PHYSICS GOALS 

ICANOE is an underground detector capable to 
achieve the full reconstruction of neutrino (and antineu- 
trino) events of any flavor, and with an energy ranging 
from the tens of MeV to the tens of GeV, for the rele- 
vant physics analyses. No other combinations can pro- 
vide such a rich spectrum of physical observations, in- 
cluding the systematic, on-line monitoring of the CNGS 
V-beam at the LNGS site. The unique lepton capabilities 
of ICANOE are really fundamental in tagging the neu- 
trino flavor. In general, the oscillation pattern of the neu- 
trinos may be complicated and involve a combination of 
Vfi — > v x , V/j — ► V e and v^, — ► V sterile transitions. In or- 
der to fully sort out the mixing matrix, unambiguous neu- 
trino flavor identification is mandatory to distinguish x's 
from v x 's and electrons from V e 's interactions. In other 
words, we stress the importance of constraining the oscil- 
lation scenarios by coupling appearance in several differ- 
ent channels and disappearance signatures. 

The sensitivity in the classic (sin' 28, Am 2 ) plot is evi- 
denced in Figure 2, for a data taking time of 4 years, with 
4.5 x 10 19 pot at each year. We remark: 

1. The recent results on atmospheric neutrinos ("A" 
and "B" of Figure 2) can be thoroughly explored 
by appearance and disappearance experiments. For 
the current central value, both CNGS and cosmic 
ray data will give independent and complemen- 




FIGURE 2. Overview of the status of the neutrino oscillations 
searches, displayed assuming two neutrino mixing schemes in 
the (sin 26, Am 2 ) plane. The 90%C.L. allowed regions ob- 
tained from the Kamiokande (resp. Superkamiokande) FC 
and PC samples are shown as A) (resp. B)). The 90% (resp. 
99%C.L.) regions consistent with the LSND excess are shown 
as dark (resp. light) shaded areas in the upper region of the 
plane. The shaded area in the region Am 2 « 10 -5 eV 2 repre- 
sents the large angle MSW solution of the solar neutrino deficit. 
CHORUS and NOMAD 90%C.L. limits on -> v t oscillations 
are visible in the upper Am 2 region. The ICANOE sensitivities 
at 90%C.L. are indicated by three curves: the limit by direct ob- 
servation of the atmospheric neutrinos (''ICANOE atm"); the di- 
rect tau appearance search at the CNGS ("ICANOE v^, — » v t "); 
the direct electron appearance search at the CNGS ("ICANOE 
-> V e "). 



tary measurements and they will provide a precise 
(sin 2 20, Am 2 ) determination. 

2. In the mass range of LSND, the sensitivity is suffi- 
cient in order to solve definitely the puzzle. 

3. At high masses of cosmological relevance for 
Am 2 <» 10 eV 2 , the sensitivity to — ► V x oscil- 
lations is better or equal to the one of CHORUS and 
NOMAD. 

4. In the atmospheric neutrino events, one can reach a 
level of sensitivity sufficient to detect also the effect 



until now observed in solar neutrinos. This purely 
terrestrial detection of the LMA solar neutrino solu- 
tion is performed using neutrino in the GeV range, 
much higher than the one of solar neutrinos. 

Since we can observe and unambiguously identify 
both v e and v x components, the full (3 x 3) mixing ma- 
trix can be explored. By itself, this is one of the main 
justifications for the choice of the detector's mass. 

In the cosmic ray channel, all specific modes (elec- 
tron, muon, NC) are equally well observed without detec- 
tor biases and down to kinematical threshold. The CR- 
spectrum being rather poorly known, a confirmation of 
the SuperKamiokande result requires detecting both (1) 
the modulation in the muon channel and (2) the lack of 
effect of the electron channel. The consistency of the 
simultaneous observation of the L/E phenomenon in as 
many modes as they are available is a powerful tool in 
separating genuine flavour oscillations from exotic sce- 
narios. 

In some favourable conditions, the direct appearance 
of the oscillated tau neutrino may be directly identified in 
the upgoing events, since even a few events will be highly 
significant. 

While in atmospheric neutrinos, the knowledge of the 
sign of the muon is of little relevance, in the case of the 
CNGS is a powerful tool to verify the neutrino nature 
after oscillation path, excluding for instance oscillation 
channels into anti-neutrinos. 

For a discussion on the nucleon decay searches, see 
Ref. (7). 



PHYSICS AT THE CNGS 

The design and performance of the CERN neutrino 
beam to Gran Sasso - the CNGS facility - are described 
in a conceptual technical design report (8). The CNGS 
beam performance for the new reference beam are sum- 
marized in Table 1 . The rms radius of the CC event 
distribution is about 1 .37 km at Gran Sasso. The expected 
numbers of detectable v x for sin 2 28= 1 and a few typical 
values of Am 2 are shown in Table 2. 

Events will occur in the whole ICANOE detector. As 
reference, we assume an exposure of 20 kton x year for 
the liquid argon. This corresponds to four years running 
of the CNGS beam in shared mode. For the events oc- 
curing in the solid detector, given the smaller mass, the 
reference exposure is 10 kton x year. The last three me- 
ters of the liquid target are defined as a transition region, 
since beam events occuring in this region are most likely 
to deposit energy in both targets. Table 3 shows the com- 
puted total event rates for each neutrino species present in 
the beam for the liquid, solid and in the transition region. 



Table 1. Predicted performance of the new CNGS reference 
beam for an isoscalar target. The statistical accuracy of the 
Monte-Carlo simulations is 1 % for the component of the 
beam, somewhat larger for the other neutrino species. 



Energy region Ey [GeV] 


1-30 


1-100 


Vfi [irT 2 /pot] 


7.1 x 10~ 9 


7.45 x 10~ 9 


Vfj CC events/pot/kt 


4.70 x 10~ 17 


5.44 x 10~ 17 


influence [ GeV ] 




17 


fraction of other events: 








o.s 


% 




2.0% 


Ve/Vf, 


0.05 % 



Table 2. Expected number of v x CC events at 
Gran Sasso per kton per year for an isoscalar 
target. Results of simulations for different val- 
ues of Am 2 and for sin 2 (28) = 1 are given for 
4.5 x 10 19 pot/year. These event numbers do 
not take detector efficiencies into account. 



Energy region E Vz [GeV] 


1-30 


1-100 


Am 2 = 1 x 10~ 3 eV 2 


2.34 


2.48 


Am 2 = 3 x 10~ 3 eV 2 


20.7 


21.4 


Am 2 = 5 x 10~ 3 eV 2 


55.9 


57.7 


Am 2 = 1 x 10~ 2 eV 2 


195 


202 



Table 3 also shows, for three different values of Am 2 , the 
V x CC rates expected in case oscillations take place. 

Event kinematics and tau identification 

Kinematical identification of the x decay, which fol- 
lows the v x CC interaction requires excellent detector per- 
formance: good calorimetric features together with track- 
ing and event topology reconstruction capabilities. The 
background from standard processes are, depending on 
the decay mode of the tau lepton considered, the v e CC 
events and/or the v^, CC and v NC events. 

In order to separate separate v x events from the back- 
ground, two basic criteria, already adopted by the short 
baseline NOMAD experiment, can be used: 

• an unbalanced total transverse momentum due to 
neutrinos produced in the x decay, 

• a kinematical isolation of hadronic prongs and miss- 
ing momentum in the transverse plane. 

In addition, given the baseline L between CERN and 
GranSasso, for the lower Aw values of the allowed re- 



Table 3. Expected event rates for an exposure of 20 kton x 
year for the liquid target and 10 ktonx year for the solid tar- 
get. All the rates include nuclear corrections and are com- 
puted for the proper target composition. For standard pro- 
cesses, no oscillations is assumed. For v x CC, we take two 
neutrino v^, — * v x with sin 2 29 = 1. 



Process 


liquid target 


transition 


solid 



v^CC 


54300 


10200 


27150 


%CC 


1090 


200 


545 


We CC 


437 


80 


219 


V e CC 


29 


5 


15 


vNC 


17750 


3330 


8875 


vNC 


410 


77 


205 



v x CC, Am 2 (eV 2 ) 









1 x 10~ 3 


52 


10 


26 


2x 10~ 3 


208 


40 


104 


3.5 x 10~ 3 


620 


115 


310 


5 x 10~ 3 


1250 


235 


625 


7.5 x 10~ 3 


2850 


535 


1425 


1 x 10~ 2 


4330 


810 


2165 



gion indicated by the atmospheric neutrino results, we ex- 
pect most of the oscillation to occur at low energy. In this 
case, a criteria on the visible energy is also very important 
to suppress backgrounds. 

In order to apply the most efficient kinematic selec- 
tion, it is mandatory to reconstruct with the best possible 
resolution the energy and the angle of the hadronic jet, 
with a particular attention to the tails of the distributions. 
Therefore, the energy flow algorithm should be designed 
with care taking into account the needs of the tau search 
analyses. 

A specially developed energy flow algorithm has been 
tested on a sample of fully simulated v e CC events, in 
order to estimate the resolution of the kinematical recon- 
struction on realistic events. It yields an average miss- 
ing Pt of 450 MeV/c. This value improves to an aver- 
age of 410 MeV/c when the primary vertex is required 
to lie within a fiducial volume of transverse dimensions 
7.8 x 7.8m 2 . 

We used the neutrino data collected in the NOMAD 
detector to probe the reliability of the physics simulation. 

CC events have been fully simulated and reconstructed 
using NOMAD official packages. We found that the kine- 
matics in the transverse plane are well reproduced by the 
Monte-Carlo model. This is clearly not the case when 
nuclear corrections are neglected. 



Table 4. Rejection of the v f CC background in the x — » e analysis. Figures are normalized to 
an exposure of 20 kton x year. 



Cuts 


v x Eff. 




V c 


v x CC 


v x CC 


v x CC 




(%) 


CC 


CC 


Am 2 = 


Am 2 = 


A/77 2 = 










10~ 3 eV 2 


3.5 x 10~ 3 eV 2 


10~ 2 eV 2 


Initial 


100 


437 


29 


9.3 


111 


779 


Fiducial volume 


88 


383 


25 


8.2 


97 


686 


One candidate with 














momentum > 1 GeV 


72 


365 


25 


6.7 


80 


561 


E vis < 18 GeV 


67 


64 


5 


6.2 


75 


522 


Pf < 0.9 GeV 


54 


31 


3 


5.0 


60 


421 


P 1 * p > 0.3 GeV 


51 


29 


2 


4.7 


56 


397 


Pf ss > 0.6 GeV 


33 


4 


0.4 


3.1 


37 


257 



— > v T appearance searches 

The channel of tau decaying into an electron plus 
two neutrinos provides the best sample for v x appearance 
studies due to the low background level. The intrinsic v,,, 
V e contaminations of the beam amount to « 470 events 
for an exposure of 20 kton x year. 

The comparison of this figure with the expected num- 
ber of v x CC events decaying into electrons shows that the 
search of x — > e at the CNGS will have to be optimized a 
posteriori. Indeed the v x rate has a strong dependence on 
the exact value of the Am 2 in the parameter region sug- 
gested by the Super-Kamiokande data, and the Am 2 value 
is not well constrained by the atmospheric neutrino ex- 
periments. 

For "large" values of Am 2 , i.e. Am 2 > 5 x 10~ 3 , the 
rate of tau is spectacular and exceeds the number of in- 
trinsic beam v e , V e CC events, i.e. S/B > 1 even prior to 
any kinematical cuts. So the kinematical cuts can be very 
mild. An excess will be striking. 

For our "best" value taken from atmospheric neutrino 
results, i.e. Am 2 = 3.5 x 10~ 3 eV 2 , the number of v x CC 
with x — > e is about 110, or about a signal over back- 
ground ratio of 1 10/470 ~ 1/4. Here with modest kine- 
matical cuts, we can extract statistically significant sig- 
nals, as shown in the following sections. 

The most difficult region lies below Am 2 « 1.5 X 
10~ 3 eV 2 , for which, kinematical cuts are tuned to sup- 
press backgrounds by a factor more than 200 while keep- 
ing about half of the signal events. 

In the following paragraphs, we discuss background 
sources and their suppression. 

V e CC rejection: The main background from genuine 
leading electrons comes from the CC interactions of the 
V e and v e components of the beam. In Table 4 we summa- 
rize the list of sequential cuts applied to reduce the V e and 



Table 5. v^ NC background to the x — > e analysis. Results are 
normalized to an exposure of 20 kton x year. We illustrate 
background reduction by means of kinematical criteria only. 
Imaging and dE jdx measurements reduce the NC background 
to a negligible level. 



Cuts 






NC 




Initial 




17750 




Fiducial volume 




15550 






Dalitz 


y conv. 


71 — > e 


7^/71° 


One candidate 


275 


4262 


6.5 


25 


P e > 1 GeV 


79 


1361 


6.3 


16 


Evh < 18 GeV 


49 


835 


3.2 


11 


Pf < 0.9 GeV 


46 


794 


1.8 


9 


P^ p > 0.3 GeV 


24 


429 


1.7 


8 


Pf ss > 0.6 GeV 


19 


350 


1.3 


7 


Imaging and dE/dx 


< 1 


< 1 


< 1 


< 1 



V e CC backgrounds and the expected number of signal 
events for three different Am 2 values. The most sensitive 
analysis predicts, for a 20 kton x year exposure, a total 
background of 4.4 events for a total x selection efficiency 
of 33%. 

V NC rejection: Neutral current events contribute to 
the background from four sources: (1) electrons from 
Dalitz decays, (2) early photon conversions, (3) interact- 
ing charged pions and (4) 7C ± /7t° overlap. Table 5 sum- 
marizes the rejection power of kinematics criteria for the 
four sources that contribute to v NC background. The re- 
quirement on the electron candidate energy E e > 1 GeV 
suppresses about one third of the Dalitz, pion overlap and 
7C° conversions induced backgrounds, since electrons in 
the jet are soft. 



20 kton x year exposure 




Missing P T Cut (GeV) 

FIGURE 3. Number of v e , V e background and tau events ex- 
pected for a 20 kton x year exposure as a function of the miss- 
ing transverse momentum. 



The ultimate discrimination of these backgrounds re- 
lies primarily on the imaging capabilities and on dE / dx 
measurements. The combination of dE /dx information 
together with kinematics criteria is sufficient to reduce 
V NC background to a negligible level. 

CC rejection: Charged current events can con- 
tribute to the background in a similar way as the neutral 
current events described above when the leading muon 
escapes detection. In case the muon is not identified, 
the event will appear in first instance as a neutral cur- 
rent event. The source of electrons which can induce 
backgrounds are then similar to those discussed previ- 
ously and are reduced to a negligible level for reasons 
already discussed. A more important source of back- 
ground specific to charged current interactions comes 
from the decays of charmed mesons. At the CNGS en- 
ergies o(v M N — > /J c X)/o(V/jN — > jtiX) « 4%, therefore 
for a total exposure of 20 kton x year we expect to col- 
lect about 200 events where a charmed meson decays into 
a positron and a neutrino. These events resemble kine- 
matically the real v x events, since they have a neutrino in 
the final state and possess a softer energy spectrum and a 
genuine sizeable missing transverse momentum. After all 
cuts, the expected number of charm induced background 
events n b cc r charm ^ for a total exposure of 20 kton x year 
is at the level of 1 event. 



Table 6. t — > e analysis summary. For a total exposure of 20 
x year we show the expected number of T events for different 
Am 2 values. The last three columns show the expected back- 
ground. 



Am 2 (eV 2 ) 


v x CC 


\ e , V e CC 


VfJ , % CC 


v^NC 


1 x 10" 3 


3 








2x 10~ 3 


12 








3x 10~ 3 


26 








3.5 x 10~ J 


35 


4.1 


1.0 


< 1 


5 x 10~ 3 


71 








7x 10~ 3 


121 








1 x 10~ 2 


248 









Combined v M — > v T sensitivity 

Table 6 summarizes the expectations for the x — > e 
analysis once kinematics criteria and muon vetoes have 
been applied to every potential background source. In 
conclusion, we obtain for a 20 kton x year exposure, that 
the overall electron selection efficiency is 32% for an ex- 
pected number of about five background events. The ex- 
pected number of fully identified tau events at the central 
Am 2 value of 3.5 x 10~ 3 eV 2 is 35. 

We show in Figure 3 as a function of the cut on the 
missing transverse momentum, the number of expected 
background and tau events for two different Am~ values. 
We see that even for a value as low as 1.5 x 10~ 3 eV 2 , a 
P™' ss cut above 0.6 GeV gives a S/B ratio in excess of 1. 

Finally it is crucial to study the exposures needed to 
obtain a statistically significant x appearance for different 
neutrino oscillation scenarios. We see in Figure 4 that for 
Am 2 = 3.5 eV 2 , few months of data taking will suffice to 
claim a 3a effect. However for Am 2 values of about 10~ 3 
eV 2 , exposures above 20 kton x year are needed to obtain 
an effect in excess of 2a. We conclude that after four 
years of running of CNGS in shared mode or after one 
year of running in dedicated mode, the ICANOE detector 
will observe a statistically significant — ► V x oscillation 
signal for most of the Am 2 values presently favored by 
atmospheric neutrino data. 



Vi, — > v e oscillation search 

The unambiguous detection and identification of v e 
CC events endows ICANOE with the ability of perform- 
ing also a — * v e oscillation search. In this case, the 
oscillation reveals itself as an excess on the number of ex- 
pected events having a leading identified electron. Given 
the expected rates for a 20 kton x year exposure, the sta- 
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FIGURE 4. x appearance statistical significance for different 
Am 2 values as a function of the exposure. The arrows indicate 
the statistical significance achieved for one year of dedicated 
run of the CNGS and four years of CNGS running in shared 
mode. 



tistical error is about 4%. Therefore the sensitivity to 
— > V e oscillations is dominated by the systematic er- 
ror on the beam knowledge. 

Figure 5 shows the 90% C.L. contours in case no os- 
cillations are observed assuming overall systematic errors 
of 5% and 10%. We observe that nearly the whole region 
favored by the LSND claim is comfortably covered. The 
excluded values are sin 2 28 > 2 x 10~ 3 at high Am 2 and 
Ain 2 > 4 x 10~ 4 eV 2 for maximal mixing. 



ATMOSPHERIC NEUTRINOS 

The physics goals of the new atmospheric neutrino 
measurements are to firmly establish the evidence of neu- 
trino oscillations with a different experimental technique, 
possibly free of systematic biases, measure the oscillation 
parameters and clarify the nature of the oscillation mech- 
anism. ICANOE will provide, in addition to comfortable 
statistics, an observation of atmospheric neutrinos of a 
very high quality. Unlike measurements obtained up to 
now in Water Cherenkov detectors, which are in prac- 
tice limited to the analysis of "single-ring" events, com- 
plicated final states with multi-pion products, occurring 
mostly at energies higher than a few GeV, will be com- 
pletely analyzed and reconstructed in ICANOE. This will 
be a significant improvement with respect to previous ob- 
servations. 
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FIGURE 5. ICANOE 90% C.L. exclusion region in case no 
V/j — * v c are experimentally observed. 



We have considered the following three methods: 

• disappearance: detection of the oscillation pat- 
tern in the L/E distribution, where L is the neutrino 
pathlength and E its energy; 

• v x appearance : comparison of the NC/CC with ex- 
pectation; 

• direct v x appearance : comparison of upward and 
downward rates of "tau-like" events. 

together with the well established ones: 

• the double ratio, (v jU /v e ) te/(v (U /v e )A/c); 

• up/down asymmetry; 

The tau appearance measurements can shed light on 
the nature of the oscillation mechanism, by discriminat- 
ing between the hypothesis of oscillations into a sterile or 
a tau neutrino. The v x appearance method is based on v x 
CC interactions with v x decaying into hadrons, hence to 
"neutral-current-like" events of high energy. An excess of 
"NC-like" events from the bottom will indicate the pres- 
ence of oscillation to the v x flavour. A kinematical anal- 
ysis of the final state particles in the event can be used to 
further improve the statistical significance of the excess. 
Such a feature can only be obtained in a detector with 
the resolution of the ICANOE liquid target, in which all 
final state particles can be identified and precisely mea- 
sured. The kinematical method would allow the evidence 
for "tau-like" events in the atmospheric neutrino beam. 

Both the Vu disappearance and the direct v x appear- 
ance methods are weakly depending on the predictions of 



Table 7. Expected atmospheric neutrino rates per ktonxyear on 
Argon in case of no oscillations. The figures have been com- 
puted with FLUKA 3D atmospheric neutrino flux with geo- 
magnetic cutoffs and include all nuclear effects. All NC events 
are included, even though only a fraction of quasi-elastic ones 
will end up in an observable final state. 



~ 80 

c 

01 

E 

.= 70 

a 
c 
o 

O 60 















50 


Process 


elastic 


single-7l 


inelastic 


Total 


<E> 














(GeV) 


40 


v„CC 


66.7 


15.9 


24.4 


107.0 


2.36 


30 


v>,CC 


12.2 


5.3 


9.8 


27.2 


3.34 




v e CC 


39.4 


8.4 


12.1 


59.9 


1.60 


20 


v e CC 


5.4 


2.1 


4.2 


11.7 


2.36 




vNC 


42.9 


8.6 


13.2 


64.8 


1.94 


10 


VNC 


21.1 
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5.0 


29.6 


2.00 
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neutrino event rates, since they rely on the comparison 
of rates induced by a downward going and upward going 
neutrinos. 

The NC/CC method, already investigated by Su- 
perKamiokande, can be significantly improved compared 
to this latter measurement. In ICANOE, imaging in the 
liquid target provides a clean bias free identification of 
neutral-current, independent on the hadronic final state, 
since the identification is based on the absence of an elec- 
tron or a muon in the final state. 

In the following sections, we will study our results for 
three different exposures: 5 ktonxyear corresponding to 
1 year of operation, 20 kton x year for 4 years and an ulti- 
mate exposure of 50 ktonxyear or 10 years of operation. 

The computed rates for the different neutrino pro- 
cesses (in events/kton/year) and their mean energies are 
quoted in Table 7, using the FLUKA-3D atmospheric 
neutrino fluxes(9). 



Event containment and muon measurement 

The muon measurement is crucial to most atmospheric 
neutrino analyses. In ICANOE, we achieve the required 
performances using the multiple scattering measurement 
rather than resorting to a high-density, coarser resolution 
detector. Keeping a low density detector, high granularity 
detector imaging allows in addition the identification and 
measurement of electrons and individual hadrons in the 
event. 

"Fully contained events" are those for which the vis- 
ible products of the neutrino interaction are completely 
contained within the detector volume. "Partially con- 
tained events" are CC events for which the muon exits 
the detector volume (only muons are penetrating enough). 



FIGURE 6. Integral distributions showing the containment for 
V^, CC as a function of the neutrino energy and the leading muon 
momentum (solid lines). Differential distribution showing the 
muon momentum resolution as a function of muon momentum 
(dashed line), including both contained and partially contained 
events. 



Figure 6 shows containment of charged current events 
for different incoming neutrino energies or muon mo- 
mentum thresholds. Clearly, because of the average en- 
ergy loss of the muon in argon (about 210 MeV/m for a 
m.i.p.), muons produced in neutrino events are often en- 
ergetic enough to escape from the subdetector volume. 

It should first be noted that for contained events the 
muon energy resolution is 4% from dE/dx measure- 
ments. For escaping muons, the high granularity of the 
imaging allows to collect a very precise determination 
of the track trajectory. Therefore the multiple scattering 
method can be effectively used to estimate the momentum 
of the escaping muons. This method requires in practice 
tracks in excess of 1 meter and works extremely well in 
the relevant energy range of atmospheric neutrino events 
(typically below 10 GeV). 

The average muon momentum resolution as a function 
of the energy threshold is shown in Figure 6. This reso- 
lution has been computed using the range measurement 
for contained muons and multiple scattering method for 
the escaping ones. For energies below 1 GeV, the average 
muon momentum resolution is about 10%. It increases 
slowly as a function of the muon momentum and reaches 
about 30% at 5 GeV. 
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FIGURE 7. Zenith angle resolution. The top plot shows the 
resolution obtained by reconstructing the incoming neutrino di- 
rection using all particles momenta, the bottom plot shows the 
resolution obtained using only leading lepton momentum. 
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Am 2 = 3.5x1 CT 3 (ev 2 ) 

± 1 a (5 Kton x year exposure, contained) 

± 1 a (50 Kton x year exposure, contained) 



Am' (eV) 

FIGURE 8. Expected R value as a function of Am 2 . Assuming a 
true A;;? 2 value of 3.5 x 10~ 3 eV 2 , the 68% confidence intervals 
are given for 5 and 50 kton x year. 



Incoming neutrino angular resolution 

The reconstruction of the zenith angle of the incom- 
ing v is of great importance in the search of oscillations 
in atmospheric neutrinos. ICANOE allows for a good re- 
construction of the incoming neutrino variables (i.e. in- 
cidence angle, energy) by using the information com- 
ing from all particles produced in the final state. Fig- 
ure 7 shows the distribution of the difference between 
the real and reconstructed neutrino angle for events with 
E v > 1 GeV. The improvement on the angular resolution 
is visible. The RMS of the distribution improves from 
~ 16 to ~ 8 degrees after the inclusion of the hadronic jet 
in the reconstruction. 



The flavor ratio and up/down asymmetry 

Given the clean event reconstruction of ICANOE, the 
ratio R of "muon-like" to "electron-like" events can be 
determined free of experimental systematic errors. In 
fact, the expected purity of the samples is above 99%. 
In particular, the contamination from n° in the "electron- 
like" sample is expected to be completely negligible. The 
measurement accuracy will be dominated by the statisti- 
cal uncertainty and by the theoretical systematic error on 
the double ratio. 

In order to estimate statistical sensitivities, we show in 
Table 8 the values and statistical errors of R for different 
exposures, assuming an oscillation — > v^, with pa- 



rameters sin 2 20 = 1 and Am 2 = 3.5 x 10~ 3 eV 2 . The ta- 
ble lists the expected results when using all events or only 
fully contained events which turn out to be quite similar. 

Clearly, after an exposure corresponding to about four 
years of running of ICANOE, the statistical error reaches 
a level below 5%. We expect that further theoretical im- 
provements should reduce the systematic error to a level 
matched to statistical precision achievable in ICANOE. 

Table 8. R as a function of the exposure assuming maximal 
mixing and Am 2 = 3.5 x 1CP 3 eV 2 . In left column: results 
obtained by using all events; right column: results obtained 
by using only fully contained events. Quoted errors are of 
statistical nature. 



Exposure (kton x year) 


R 

all events contained 


5 


0.696 ±0.067 


0.674 ±0.074 


20 


0.696 ±0.033 


0.674 ±0.037 


50 


0.696 ±0.021 


0.674 ±0.023 



From the value of R and from its zenith angle depen- 
dence we can obtain the allowed parameter regions of 
neutrino oscillations. Figure 8 shows how precisely we 
can determine Am 2 in the oscillation case. The la re- 
gions corresponding to 5 and 50 kton x year have been 
computed using contained events only. 

The expected up/down asymmetries are shown in Fig- 
ure 9 for three different mixing angles as a function of 
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FIGURE 9. Expected Up-Down asymmetry for sin 2 (29) = 1,0.8,0.6 as a function of Am 2 for low momentum (P<0.4 GeV) 
events, (left) when all particles are used to reconstruct the incoming neutrino direction (right) only the lepton is used. The error 
bands show the 1 a uncertainty for 20 ktonxyear and 50 ktonxyear exposures and a 5% systematic error assuming sin 2 26 = 1 and 
Am 2 = 3.5 x 10~ 3 eV 2 . 



Am . Note in Figure 9 the much better asymmetry reso- 
lution of ICANOE for low energy muons when compared 
to a measurement including the lepton only. 



v A , disappearance -L/E studies 

In order to verify that atmospheric neutrino disappear- 
ance is really due to neutrino oscillations, an effective 
method consists in observing the modulation given by the 
characteristic oscillation probability: 

p (Jj^j = 1 -sin 2 (29) sin 2 ^1.27Am 2 ^J (1) 

with L in km, E in GeV, Am 2 in eV 2 . This modulation 
will be characteristic of a given Am 2 , when the event rate 
is plotted as a function of the reconstructed L/E of the 
events when compared to theoretical predictions. The 
ratio of the observed and predicted spectra has the ad- 
vantage of being quite insensitive to the precise knowl- 
edge of the atmospheric neutrino flux, since the oscilla- 
tion pattern is found by dips in the L/E distribution while 
the neutrino interaction spectrum is known to be a slowly 
varying function of L/E. Such a method is in principle 
capable of measuring Am 2 exploiting atmospheric neu- 
trino events. 



A smearing of the modulation is introduced by the fi- 
nite L/E resolution of the detection method. Precise mea- 
surements of energy and direction of both the muon and 
hadrons are therefore needed in order to reconstruct pre- 
cisely the neutrino L/E, This is quite well achieved in 
ICANOE. The contained muons can be measured with 
a resolution of 4%, while the non-contained muons are 
measured by multiple scattering method . 

The RMS reconstructed L/E resolution is about 30% 
for events with E V i s n,i e > 1 GeV. 

The v A , survival probability as a function of L/E let us 
determine the value of Am~ in case of oscillation is con- 
firmed. In figure 10 we can see the survival probabilities 
of V/j for neutrino oscillation hypothesis and four differ- 
ent values of Am 2 . The first minimum on the survival 
probability happens at highest L/E values for the lowest 
Am 2 values, and allows us to discriminate between them 
for an exposure of 50 ktonxyear. 

The most favored solution of the atmospheric neutrino 
anomaly is through — > V x oscillations. However, al- 
ternative explanations, like neutrino decay, cannot yet be 
excluded (10). For example, in a model in which one of 
the mass-eigenstates of neutrinos with flavour content 
decays, the disappearance probability can be described by 
the expression: 

P(Vfi -> V. v ^) = (sin 2 9 + cos 2 0e- aL / 2£ ) 2 . (2) 




2.5 3 3.5 4 
log(L/E (km/GeV)) 



2.5 3 3.5 4 
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FIGURE 10. Survival probability as a function of the L/E ratio 
assuming neutrino oscillation hypothesis and for various Am 2 
values and for 50 ktonxyear. Only statistical error has been 
considered. 



Such a model gives an equally good fit for the choice of 
parameters a = 1/63 GeV/km, cos 2 9 = 0.30 (10). 

The capability of distinguishing between the two hy- 
pothesis depends on the resolution in measuring the L/E 
ratio, which depends on the angular and momentum res- 
olution. 

Figure 1 1 shows the survival probabilities as a func- 
tion of L/E for the neutrino decay hypothesis with a = 
ttiy/Zy = 1/63 GeV/km and cos 2 9 = 0.30, and the oscil- 
lation hypothesis with sin" 29= land Am 2 = 1.0 x 10~ 3 , 
for an exposure of 50 ktonxyear. Both hypothesis are 
distinguishable from each other at around 2000 km/GeV 
within the statistical errors. 



(Direct) appearance of tau neutrinos 

To discriminate between — > V x and — > v s oscil- 

lations, we measure the ratio R^c/e = nc cx p /v' cc"p • ^ n 
oscillation to an active neutrino leads to R^c/e — 1> while 
R^c/e ~ 0.7 is expected for an oscillation to sterile neu- 
trino. 

Table 9 shows values and errors of R^c/e m case °f 
oscillation to a sterile neutrino, for all events and fully 
contained events respectively. 

For Am 2 < 10 -2 eV 2 , oscillations of into v x would 
in fact result in an excess of "neutral-current-like" events 
produced by upward neutrinos with respect to downward, 
since charged-current v x interactions would contribute to 
the "neutral-current-like" event sample, due to the large x 



FIGURE 11. Survival probability as a function of the L/E ratio 
for oscillation (triangles) and decay hypothesis (circles). Only 
statistical error has been considered. 



Table 9. R^c/e as a function of the exposure assuming os- 
cillation to a sterile neutrino. Quoted errors are of statistical 
nature. 



Exposure (ktonxyear) 


R NC/e 

all events contained 


5 


0.674 ±0.086 


0.670 ±0.087 


20 


0.674 ±0.043 


0.670 ±0.044 


50 


0.674 ±0.027 


0.670 ±0.028 



branching ratio into hadronic channels. Moreover, due to 
threshold effect on x production, this excess would be im- 
portant at high energy. Oscillations into a sterile neutrino 
would instead result in a depletion of upward muon-less 
events. Discrimination between — > V x and — > \ s is 
thus obtained from a study of the asymmetry of upward to 
downward muon-less events. Because this method works 
with the high energy component of atmospheric neutri- 
nos, it becomes effective for relatively large values of 
Am 2 (> 3 x 10~ 3 eV 2 ). 

Charged current V x rates for five Am 2 hypothesis: 
5 x 10~ 4 , 1 x 10" 3 , 3.5 x 10~ 3 eV 2 , 5 x 10~ 3 eV 2 and 
1 x 10" 2 eV 2 are listed in Table 10. We see that the rates 
saturate at about one event per ktonxyear for the larger 
Am 2 values. Such small rates pose a major experimental 
challenge in the detection of v x in the cosmic ray induced 
neutrino flux. 

The total visible energy (E v i s u,i e ) is a suitable discrim- 
inant variable to enhance the S/B ratio. After cuts, sur- 
viving events are classified as: rij, (number of expected 



Table 10. Expected v x , V t absolute rates for five different Am 2 
with FLUKA-3D fluxes and relative to FLUKA-1D and Bartol 
fluxes. 



V T + V x CC (NUX, Fluka 3D flux) 


Rel. to 
Fluka ID 


Rel. to 
Bartol 


Am 2 (eV 2 ) 


Rate (ktonx year) 


DIS 


QE 


Sum 


5 x 1CT 4 


0.11 


0.11 


0.22 


0.96 


0.81 


1 x 1(T 3 


0.28 


0.18 


0.46 


1.02 


0.84 


3.5x 10~ 3 


0.59 


0.21 


0.80 


1.00 


0.81 


5 x 1(T 3 


0.64 


0.24 


0.88 


1.01 


0.80 


1 x 10~ 2 


0.70 


0.20 


0.90 


0.99 


0.78 



downward going background) and n t = rib + n s (number 
of expected upward going events, where n s is the number 
of taus). The statistical significance of the expected n s 
excess is evaluated following two procedures: 

• The ft, and /, pdf 's are integrated over the whole 
spectrum of possible measured r values and the 
overlap between the two is computed: P a = 
JJ° min(fb(r), f t (r))dr, where ft, and/ r are the Pois- 
son p.d.f.'s for means /j = rib and /j = «, respectively. 
The smaller the overlap integrated probability (P a ) 
the larger the significance of the expected excess. 

oo e~"bn r 

• computing the probability = J n — jr^dr that, 
due to a statistical fluctuation of the unoscillated 
data, we measure n t events or more when rib are ex- 
pected. 

For a 50 ktonx year exposure, the results of a search 
based on E^Me are shown in Table 1 1 . We see that a cut 
on visible energy between 6 and 7 GeV results in: (1) an 
overlap integrated probability between the two distribu- 
tions amounting to 25 — 26%. (2) a Poisson probability 
that the measured excess ("x bottom") corresponds to a 
statistical fluctuation is 0.6 — 0.8%. 

The search for v x appearance can be improved taking 
advantage of the special characteristics of v x CC and the 
subsequent decay of the produced x lepton when com- 
pared to CC and NC interactions of and v^, i.e. by 
making use of Pi epton and Phadron- 

The information related to the directionality of the in- 
coming neutrino (i.e. the beam direction!) is missing. 
As a result, we have three kinematical independent vari- 
ables in order to separate signal from background. After 
a careful evaluation of the performance of different com- 
binations of variables, we decided to use: E V j S n,i e , yt,j (the 
ratio between the total hadronic energy and E V i S ibi e ), and 
Qt (the transverse momentum of the x candidate with re- 
spect to the total measured momentum) which contains 



Table 11. Number of NC and tau events as a function of the 
visible energy cut. The statistical sample used corresponds 
to an exposure of 50 ktonx year. 



50 kton x year exposure 



Evisible cut 


v NC top 


x bottom 


Pa (%) 


Pf, (%) 


> 1 GeV 


327 


22 


55.0 


10.8 


> 2 GeV 


150 


22 


38.6 


3.54 


>3GeV 


95 


21 


30.6 


1.6 


>4GeV 


67 


20 


25.3 


0.8 


>5GeV 


51 


17 


27.3 


0.9 


> 6 GeV 


40 


16 


24.6 


0.6 


>7GeV 


33 


14 


26.6 


0.8 


> 8 GeV 


28 


13 


26.7 


0.8 


> 9 GeV 


23 


12 


26.2 


0.7 


> 10 GeV 


21 


11 


28.3 


0.9 



1 NC 

1 .... 1 , . .": - i-7'.. 


T 


3-2-101 2 3 


4 5 6 


r * . _j-,..^r~T . ... i .... i . . ."."i---' r - r .. 


T 


4-3-2-1 1 


2 3 4 




"1 T 



-4 -3-2-10 1 2 

FIGURE 12. Likelihood ratio distributions, arbitrarily normal- 
ized, for signal and background events. The distributions are 
computed separately for each of the considered hadronic chan- 
nels: 71, p, 3 prongs. 



the information on the isolation of the tau candidate from 
the recoiling jet. 

The chosen variables are not independent one from an- 
other but show correlations between them. These corre- 
lations can be exploited to reduce the background. In or- 
der to maximize the separation between signal and back- 
ground, we use three dimensional likelihood functions 
-L{QT,E v j S jbi e ,ybj) where correlations are taken into ac- 
count (see Figure 12). 

Table 12 illustrates the statistical significance achieved 



Table 12. Expected background and signal events for dif- 
ferent combinations of the it, p and 3k analyses. The con- 
sidered statistical sample corresponds to an exposure of 50 
ktonxyear. 



K 


P 


3ji 


Top 


Bot. 


Pa 




Cut 


Cut 


Cut 


Evts 


Evts 


(%) 


(%) 





0.5 





112 


134 


32.1 


1.9 (2.3a) 


1.5 


1.5 





46 


63 


24.8 


0.7 (2.7a) 


3 


-1 





43 


59 


26.1 


0.8 (2.6a) 


3 


0.5 





12 


23 


18.3 


0.14 (3.3a) 


3 


1.5 





10 


20 


18.8 


0.16 (3.2a) 


3 


0.5 


-1 


30 


45 


21.9 


0.4 (2.9a) 


3 


0.5 


1 


9 


17 


25.9 


0.5 (2.8a) 



by several selected combinations of the likelihood ra- 
tios for an exposure equivalent to 50 ktonxyear. We 
take as the best combination the one with the lowest P a . 
This is achieved for the following set of cuts: lnX K > 3, 
lnA,p > 0.5 and lnA^ > 0. The expected number of NC 
background events amounts to 12 (top) while 12+11 = 
23 (bottom) are expected. This corresponds to a P a of 
18.3%. In the case we consider E v j S jbi e as the unique 
discriminating variable, a similar number of background 
events is obtained demanding E v j s ibi e > 14 GeV. With this 
cut, the expected number of x events is 7 and the P a is 
37%. Therefore, for the same level of background, the 
approach using the ratio of three dimensional likelihood 
functions enhances the number of expected signal events 
by approximately 50%. 

Finally, in figure 13 we present the Poisson probabil- 
ity Pp for the measured excess of upward going events to 
be due to a statistical fluctuation as a function of the ex- 
posure. The bottom curve corresponds to the case where 
no kinematical selection has been applied and only a cut 
on Edible > 6 GeV is used. We see that for exposures 
around 30 ktonxyear, in case we use the kinematical se- 
lection algorithm, the observed excess corresponds to a 
2.6a effect. This effect is larger than 3o for an exposure 
of 50 kton x year. 
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FIGURE 13. Probability for the measured excess of upward 
going events to be due to a statistical fluctuation of the data as a 
function of the exposure. 
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